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Direct and phonon-assisted tunneling currents in InAlGaAs-InGaAs bulk and double quantum 
well interband tunnel heteroj unctions are simulated rigorously using the non-equilibrium Green's 
function formalism for coherent and dissipative quantum transport in combination with a simple 
two-band tight-binding model for the electronic structure. A realistic band profile and associated 
built-in electrostatic field is obtained via self-consistent coupling of the transport formalism to 
Poisson's equation. The model reproduces experimentally observed features in the current-voltage 
characteristics of the device, such as the structure appearing in the negative differential resistance 
regime due to quantization of emitter states. Local maps of density of states and current spec- 
trum reveal the impact of quasi-bound states, electric fields and electron-phonon scattering on the 
interband tunneling current. In this way, resonances appearing in the current through the double 
quantum well structure in the negative differential resistance regime can be related to the alignment 
of subbands in the coupled quantum wells. 
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Interband tunnel junctions (TJ) are essential components 
in multi-junction solar cell devices, allowing for the se- 
ries connection of single junction subcells with different 
band gaps as used for an enhanced utilization of the solar 
spectrumP^l. In a recent publication, the use of a double 
quantum well (DQW) structure in the junction region 
was shown to enhance considerably the peak tunneling 
current in an InAlGaAs-InGaAs TJ, at moderate loss of 
transparency. For a proper analysis of the potential 
of this approach and to enable further optimization of 
the junction design, an accurate and comprehensive the- 
oretical description of the physical processes involved in 
the transport is required. However, conventional models 
for tunnel junction as used for the simulation of multi- 
junction solar cells are usually limited in this respect by 
a large number of simplifying assumptions on the density 
of states participating in the tunneling, the occupation of 
these states and the assisting scattering mechanisms^. 

In this paper, the tunnel junction is simulated using a 
rigorous quantum transport approach based on the non- 
equilibrium Green's function (NEGF) formalism as used 
in the nano-electronics community for the simulation of 
interband tunneling transistors^"^. In combination with 
atomistic electronic structure theory and self-consistent 
coupling to a Poisson solver, the approach is able to pro- 
vide a realistic picture of the relevant density of states 
and associated current flow under non-equilibrium condi- 
tions, by accurately reflecting the effects of internal fields, 
carrier confinement and elastic as well as inelastic scat- 
tering mechanisms. 

The present approach is based on the implementation 
of the NEGF formalism for III-V quantum well (QW) 
structures with photvoltaic applications as introduced in 
Ref. [11 , with the material parameters of the nearest- 
neighbor two-band tight-binding (TB) model and in- 
teractions adjusted to the InAlGaAs-InGaAs system of 



Ref. [3 , as presented in Tab. [TJ The orbital basis, in dif- 
ference to a band basis like k-p, enables a consistent and 
unified description of electrons independent of their po- 
sition in energy with respect to conduction and valence 
bands and is thus ideally suited for interband tunneling 
processes. The TB parameters are adjusted to repro- 
duce the zone center band gaps and effective masses of 
the bulk component materials, choosing the light hole 
mass for the description of the valence band. Among the 
various existing scattering mechanisms for electrons and 
holes, elastic coupling to acoustic phonons is considered 
within the deformation potential approximation, and in- 
elastic interaction with polar optical phonons is included 
via the corresponding Frohlich Hamiltonian, with param- 
eters also displayed in Tab. [I[ While the former process 
results in a broadening of the spectral quantities in both 
energy and transverse momentum domain, the latter is 
essential for the description of the phonon mediated tun- 
neling current at vanishing overlap of electron and hole 
states in the injector regions. At the present stage, no 
trap-assisted tunneling relying on additional defect states 
is included, which potentially leads to an underestima- 



TABLE I. Material parameters used in simulations 





Ino.52Alo.33Gao.15As 


Ino.53Gao.47As 


E s [eV] 


1.04 


0.74 


E p [eV] 


-0.14 


0.0 


V sp [eV] 


2.54 


2.6 


mli /m 


0.065 


0.041 


mt h /mo 


0.087 


0.052 


e r 


13.2 


14.2 


Coo 


10.7 


11.0 


hVL LO [eV] 


0.038 


0.033 


D ac [eV] 
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5 
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FIG. 1. Schematic layer and material structure of the bulk 
and DQW tunnel- heteroj unction devices from Ref. [3[. The 
doping density amounts to N d / a = 2 x 10 17 
lightly doped n + /_p + regions and N d / a — 10 19 
heavy n ++ /_p ++ doping. 
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tion of the dark interband current, especially at higher 
forward bias voltages. 

The device structures under consideration are dis- 
played in Fig. [I] and correspond to the central re- 
gion of the devices discussed in Ref. [3]. The tun- 
nel junction consists of two adjacent 15 nm thick lay- 
ers of Ino.52Alo.33Gao.15As with high doping level of 
Nd/a = 10 19 cm -3 . In the DQW structure, 8 nm of 
Ino.53Gao.47As are inserted on each side of the junction 
to form the potential wells leading to a local constriction 
of the interband tunneling barrier. For the present sim- 
ulation, the size of the lightly doped (N d / a = 2 x 10 1T 
cm -3 ) buffer layers on top and bottom of the tunnel junc- 
tion is reduced to 30 nm, and neither the p-doped InGaAs 
capping layer nor the n-doped InP substrate present in 
the original structures are included. 

The equilibrium band edge profiles and associated local 
density of states (LDOS) at vanishing transverse momen- 
tum (k|| = 0) corresponding to the above layer structures 
and doping levels are shown in Fig. [2] In the bulk het- 
erojunction case, displayed in Fig. [2|a), the strong band 
bending close to the junction results in the appearance of 
bound states in both electron and hole emitter regions. 
While these states feed the tunneling current, they are 
themselves not accessible from the contacts without en- 
ergy relaxation. Hence, there is only a vanishingly narrow 
forward bias regime with finite energetic overlap of bulk 
electrode states at the low doping levels of the n + and 
p+ regions, where ballistic transport from one contact to 
the other via exclusively direct tunneling in the junction 
region is possible, in spite of the lower band gap at the 
n-contact. As a consequence, any description based on 
a ballistic picture of direct tunneling will severely under- 
estimate current flow even close to zero bias. There is, 
however, an extended forward bias regime with finite en- 
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FIG. 2. Band profile and local density of states at vanish- 
ing transverse momentum and zero bias voltage, for (a) the 
bulk heteroj unction device, exhibiting state quantization in 
the highly doped emitter regions adjacent to the junction 
and (b) the DQW heteroj unction device, where the strong 
band bending leads to the unbinding of the high-lying con- 
fined states. 



ergetic overlap of interface states close to the junction 
via which the tunneling can proceed. In the DQW case, 
displayed in Fig.[2|b), both the energetic overlap and spa- 
tial extension of interface states is increased due to the 
presence of the additional confined subgap states in the 
quasi two-dimensional QW regions close to the junction. 
The extremely strong band bending results in a partial 
unbinding of the higher QW states, such that the LDOS 
is no longer well described by the states of a square well 
potential, but rather a new confinement potential reflect- 
ing the effects of both band bending and band offsets. In 
combination with the confined states at the heteroj unc- 
tion interface, the injector regions thus provide a complex 
combination of bound, quasibound and continuum states 
participating in the interband tunneling process. 

Fig. [3] displays the current-voltage characteristics of 
the interband charge transport in the bulk and DQW 
heteroj unctions as computed with the NEGF model. The 
current-voltage characteristics reproduce the qualitative 
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FIG. 3. Current- volt age characteristics of bulk and DQW het- 
eroj unction devices at forward bias. Due to the lower effective 
barrier potential, the peak current in the DQW device exceeds 
that of the bulk junction by more than two orders of magni- 
tude and occurs at much larger bias voltage. The structure 
in the negative-differential resistance regime exhibited by the 
bulk device can be ascribed to the quasi- 2D state quantization 
in the emitter regions adjacent to the junction, as shown in 
Fig. [2] The corresponding features in the IV-characteristics 
of the DQW heteroj unction device are considerably stronger 
and are identified as the signatures of resonant tunneling be- 
tween QW states. 
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features of the experimental data in Ref. [3] in terms of 
significantly higher peak current and peak voltage, result- 
ing from the lower effective barrier potential, but lower 
negative differential resistance (NDR) for the DQW TJ 
as compared to the bulk TJ. As can be inferred from 
Fig. [2] (a), the structure appearing in the NDR region of 
the IV characteristics for the bulk heteroj unction is due 
to discretization of the longitudinal DOS in the injection 
regions adjacent to the junction. In the DQW sample, 
the local current maxima are much more pronounced and 
are identified as the signatures of resonant tunneling be- 
tween different QW states. 

In order to analyze the current flow in the two struc- 
tures in more detail with respect to the different extended 
and localized states involved, the current spectrum at dif- 
ferent bias voltage is computed for the full spatial extent 
of the device. In the case of the bulk heterojunction, the 
result displayed in Fig. [4| together with the correspond- 
ing LDOS (k|| = 0), allows the distinction of different 
regimes of interband charge transport depending on the 
applied bias: (a) At the peak voltage of 0.06 V, current 
flow is dominated by carriers first tunneling through the 
triangular barrier of the hetero-interface at the electron 
contact into the lowest subband of the doping-induced 
potential well at the n-side of the junction and then into 
the second subband of the doping-induced potential well 
at the p-side, from where escape to the hole contact pro- 
ceeds via electron-phonon scattering. A minor contribu- 
tion is due ballistic transport above the triangular bar- 
rier into the subband associated with the lowest localized 



FIG. 4. Current- voltage characteristics, local density of states 
and corresponding local current spectrum of the bulk tunnel- 
hetero junction device at forward bias voltages of (a) 0.06 V, 
(b) 0.16 V and (c) 0.26 V, illustrating different interband 
tunneling processes depending on the energetic alignment of 
localized and extended states in the injection regions adjacent 
to the junction. 



state at the p-side; (b) At a bias of 0.16 V, after initial re- 
laxation close the contacts mediated by multiple scatter- 
ing events with phonons, current flow proceeds again pre- 
dominantly via tunneling through the triangular barrier 
from the n-side, but now into the lowest subband of the 
doping-induced potential well at the p-side; (c) At even 
larger bias voltage of 0.26 V, the major part of the current 
flows via the lowest subband of the triangular potential 
well at the hetero-interface into the lowest subband at the 
p-side, overcoming a substantially thicker tunneling bar- 
rier, while only a small fraction of carriers are transported 
via phonon-assisted tunneling from the higher interface 
states. Beyond this transport regime, i.e., for consider- 
ably larger separations of contact Fermi levels, there are 
no more available current paths in terms of energetically 
overlapping states in the injection regions, and current 
breaks down exponentially, until thermionic emission sets 
in, or, in more realistic cases, defect-assisted tunneling 
takes over. In the DQW TJ, as shown in Fig. [5] the dom- 
inant role in the transport process is played by the QW 
states, which in case of energetic alignement at specific 
values of the bias voltage give rise to resonant tunnel- 
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FIG. 5. Current-voltage characteristics, local density of states 
and corresponding local current spectrum of the DQW tunnel- 
hetero junction device at forward bias voltages of (a) 0.16 V, 
(b) 0.21 V and (c) 0.31 V, at which electrons tunnel resonantly 
from the lowest QW subband at the n-side of the junction 
to the third, second and first QW subband at the p-side, 
respectively. 



ing. On the n-side, current flow proceeds over the whole 
bias regime via tunneling into the second subband of the 
QW and subsequent relaxation into the lowest subband, 



while the states in the triangular well at the heterointer- 
face seem to be irrelevant for transport. On the p-side, 
there is a strong bias dependence of the specific subband 
involved in the interband tunneling process: (a) at 0.16 
V, the alignement is between the lowest subband at the 
n-side and the third subband state at the p-side, (b) at 
0.21 V, the second subband is resonant and (c) at 0.31 
V, carriers are transferred between the lowest subbands 
of the QW. Again, beyond this point, energetic overlap 
is limited to phonon sidebands, which results in an expo- 
nential decay of the current with increasing forward bias 
voltage. Finally, in both tunnel junction devices, due to 
the prominent role of narrow bands with enhanced trans- 
mission, the current spectrum close to the contacts shows 
a pronounced structure depending on the phonon energy, 
in spite of the bulk-like, continuous DOS in these regions. 

In conclusion, we have presented a rigorous approach 
to the simulation of carrier transport in interband tunnel 
junctions with applications in multi-junction solar cells, 
which is able to account for both direct and phonon- 
assisted contributions under consideration of realistic 
spatial potential variations, including arbitrarily shaped 
quantum well structures. By consideration of the local 
maps of DOS and current spectrum, a microscopic pic- 
ture of the current flow at different values of forward 
bias voltage is obtained. Apart from providing insight 
into the dominant transport mechanisms at given bias, 
the simulation results are in qualitative agreement with 
the experimentally found large increase in peak current 
and peak voltage as well as a decrease of NDR by the 
insertion of a double quantum well structure in the junc- 
tion region. The quantitative agreement of the tunneling 
currents with the measured values is expected to be im- 
proved by the use of a larger tight-binding basis such 
as sp 3 d 5 s* 12 , a more accurate treatment of the scatter- 
ing processes and, most notably, the inclusion of defect- 
assisted tunneling processes. The latter is essential to ex- 
tend the present approach to multi-junction solar cell de- 
vices based on disordered and defect-rich materials, such 
as the micromorph silicon thin-film tandem solar cell. 
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